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first two features) and moderate to high thrust (limited primarily by the 
average laser power available)* The present research addresses questions 
related to thruster performance and optical design. j 

In the thruster scheme under consideration, parabolic nozsle mils focus 
the incoming laser beam to yield breakdown in a propellant at the focal 
point of the parabola. The resulting high pressure plasma is characteris¬ 
tic of a detonation wave initiation by high power laser-induced breakdown. 
With a short laser pulse, the detonation wave quickly becomes a blast wave 
which propagates to the nozzle exit plane converting the high pressure of 
the gas behind it to a force on the nozzle wall. Propellant is fed to the 
focal region from a plenum chamber. The laser-induced blast wave stops 
the propellant flow through the throat until the pressure at the throat 
decays to the sonic pressure; then the propellant flow restarts. The 
process is repeated with each successive laser pulse. | 

j 

A major part of the present work Involved c;trrying out thruster perfor¬ 
mance experiments using a small scale conical nozzle, external focusing 
optics, and pulsed CO 2 (X*10.6 wm) and pulsed XeF (X>0.35 imi) lasers as 
the propellant heating sources. All experiments were performed in a 
vacuum environment. With the XeF laser (Ep~2.7 J, Tp^O.d vs) and argon 
propellant, a maximum specific impulse of 500 seconds was achieved with a 
conversion efficiency of laser energy to propellant kinetic energy upwards 
of 50%. With the same propellant and the CO 2 TEA laser (Ep*‘6.5 <7, ‘tp*3 y 
ys), the maximum specific impulse was 700 seconds with a laser to kinetic 
energy conversion efficiency " 20%. Finally, using hydrogen propellant 
and the CO 2 laser as the heating source, the experimental data indicate 
specific Impulse to 3000 seconds with corresponding energy conversion 
efficiencies upwards of ■■ 70%. 

. . i 

* f 

Zn.fiupport of the truster performance experiments, a detailed computer 
aK>^l that was developed previously was modified and used to simulate the 
mo^Ic fluid dynamiqs in the conical thruster. The model treats non- 
stj(^dy nozzle flow of a real gas with laser heat addition. Major enphasis 
in the present work tias directed toward simulating the fluid dynamics ap¬ 
propriate for repetitively-pulsed operation, i.e., the situation in which 
the nozzle has only partly refilled with fresh gas between successive 
laser pulses emd the remaining downstream section of the nozzle contains 
the low pressure tail of the previous blast wave. For a given initial 
propellant mass distribution 2 md laser input, the computer program calcu¬ 
lates such quantities as the exhaust gas momentum, kinetic energy, and 
specific impulse. 

Finally, an assessment has been made of optical design criteria for a 
self-focusing paraboloidal nozzle. The primary conclusions reached ares 
(1) the mission requirements for off-axis performance allow the surface to 
be of only very modest optical quality, (2) for off-axis angles t 10°, 
only very modest magnifications 25} can be anticipated, (3) optimum op¬ 
tical and fluid dynamic performance appear to be mutually inclusive, con¬ 
sequently, a compromise is required, and (4) the impact of off-axis focus¬ 
ing and breakdown on thruster performance should be evaluated. 












ABSTRACT 




In this report we present results of continuing analytical and experi¬ 
mental Investigations carried out to evaluate the concept of pulsed laser pro¬ 
pulsion. This advanced propulsion scheme, which has been the subject of sev¬ 
eral previous studies. Involves supplying propellant energy by beaming short, 
repetitive laser pulses to a thruster from a remote laser power station. The 
concept offers the advantages of a remote power source, high specific isqpulse, 
high payload to total mass ratio {a consequence of the first two features) and 
moderate to high thrust (limited primarily by the average laser power avail¬ 
able). The present research addresses questions related to thruster perfor¬ 
mance and optical design. 

In the thruster scheme under consideration, paralsollc nozzle walls 
focus the Incoming laser beam to yield breakdown In a propellant at the focal 
point of the parabola. The resulting high pressure plasma Is characteristic 
of a detonation wave Initiation by high power laser-induced breakdown. With a 
short laser pulse, the detonation wave quickly becomes a blast wave which 
propagates to the nozzle exit plane converting the high pressure of the gas 
behind It to a force on the nozzle wall. Propellant Is fed to the focal 
region from a plenum chamber. The laser-induced blast wave stops the propel¬ 
lant flow through the throat until the pressure at the throat decays to the 
sonic pressure; then the propellant flow restarts. The process Is repeated 
%d.th each successive laser pulse. _ 

A major part of the present work Involved carrying out thruster per¬ 
formance experiments using a small scale conical nozzle, external focusing op¬ 
tics, and pulsed CO 2 (1>10.6 urn) and pulsed XeF (X>0.35 urn) lasers as the pro¬ 
pellant heating sources. All experiments were performed In a vacuum environ¬ 
ment. With the XeF laser (Ep>>2.7 J, Tp>>0.6 vs) and argon propellant, a maxi¬ 
mum specific Impulse of 500 seconds was achieved with a conversion efficiency 
of laser energy to propellant kinetic energy upwards of 50%. With the same 
propellant and the CO 2 TEA laser (Ep*6.5 J, ‘Tp''3 vs), the maximum specific 
Impulse was 700 seconds with a laser to kinetic energy conversion efficiency 
■ 20%. Finally, using hydrogen propellant and the OO 2 laser as the heating 


a 








source, the experimental data indicate specific impulse up to 3000 seconds 
with corresponding energy conversion efficiencies upwards of " 70%. 

In support of the thruster performance experiments, a detailed comput¬ 
er model that was developed previously was modified and used to simulate the 
nozzle fluid dynamics in the conical thruster. Itie model treats non-steady 
nozzle flow of a real gas with laser heat addition. Major emphasis in the 
present work was directed toward simulating the fluid dynamics appropriate for 
repetitively-pulsed operation, l.e., the situation in which the nozzle has 
only partly refilled with fresh gas between successive laser pulses and the 
remaining downstream section of the nozzle contains the low pressure tail of 
the previous blast wave. For a given initial propellant mass distribution and 
laser input, the computer program calculates such quantities as the exhaust 
gas momentum, kinetic energy, and specific impulse. 

Finally, an assessment has been made of optical design criteria for a 
self-focusing paraboloidal nozzle. The primary conclusions reached are: 

(1) the mission requirements for off-axis performance allow the surface to be 
of only very modest optical quality, (2) for off-axis angles i 10^, only very 
modest magnifications 25) can be anticipated, (3) optimum optical and fluid 
dynamic performance appear to be mutually inclusive, consequently, a compro¬ 
mise is required, and (4) the impact of off-axis focusing and breakdown on 
thruster performance should be evaluated. 
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INTRODtXTTION AND EXECUTIVE SUMMARY 


Repetitively-pulsed (RP) laser propulsion Is a pulsed jet propulsion 
concept In which the propellant energy Is supplied by the absorption of short, 
repetitive laser pulses beamed to the thruster from a remote laser power sta¬ 
tion. nils advanced propulsion concept offers the advantages of a remote 
power source, high specific Impulse, high payload to total mass ratio (a con¬ 
sequence of the first two features), and moderate to high thrust (limited pri¬ 
marily by the average laser power available). 

In the RP thruster concept shown In Fig. 1.1, parabolic nozzle walls 
focus the Incoming beam to yield breakdown In a propellant at the focal point 
of the parabola. The resulting high pressure plasma Is characteristic of det¬ 
onation wave intlatlon by high power laser-induced breakdown. With a short 
laser pulse, the detonation wave quickly becomes a blast wave, which propa¬ 
gates to the nozzle exit plane converting the high pressure of the propellant 
gas behind It into a force on the nozzle wall. Propellant Is fed to the focal 
region from a plenum chamber. The laser-induced blast wave stops the propel¬ 
lant flow through the throat until the pressure at the tluroat decays to the 
sonic pressure; then, the propellant flow restarts. The process Is repeated 
with each successive laser pulse. 

Physical Sciences Inc. (PSI) has been developing the technology of 
pulsed laser propulsion under DARPA support.^Previous experiments and 
supporting theoretical studies have shown that repetitively pulsed 10.6 urn 
laser radiation can be utilized to heat simple propellants to high stagnation 
temperatures; the result Is rocket propulsion at high specific impulse (Igp). 
Using a simple thruster concept (similar to the one shown In Fig. 1.1) with a 
variety of propellants, we obtained specific Impulses of up to 1000 s with an 
energy conversion efficiency of 40%. PSI has since been engaged In utilizing 
the knowledge and experience obtained from this program to advance laser pro¬ 
pulsion technology for DARPA space applications. Since laser propulsion ac¬ 
tivities are driven by advances In laser technology for other applications In 
space. It has been proposed that pulsed laser propulsion systems utilize vls- 
Ible/UV pulsed lasers. Therefore, the recent program objectives have been to 
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1) establish missions and associated laser requirements for a ground based 
laser propulsion system, 2) establish a thruster design capability for effi¬ 
cient performance as a function of desired specific impulse and laser wave¬ 
length, 3) develop scaling laws for laser propulsion with pulsed visible 
lasers, and 4) Identify the technology Issues for laser propulsion with vis¬ 
ible lasers. 

The technical approach for the previous program phaseInvolved per¬ 
forming theoretical modeling and systems analysis along with supporting exper¬ 
imental studies using a commercially built pulsed XeF laser. Although It has 
not been established that XeF (0.35 urn) Is Indeed the desirable wavelength, 
the advanced nature of exclmer laser technology has led to the availability of 
commercially built lasers with sufficient energy and the proper pulse time to 
perform meaningful scaling experiments. The analytical aspects of this pro¬ 
gram were divided into studies of propulsion physics wavelength scaling and 
mission analysis to identify DOD applications of laser propulsion systems. 

The propulsion physics studies were carried out at PSI, and the mission analy¬ 
sis was performed both at PSI and the Lockheed Missiles and ^ace Company 
under subcontract to PSI. The experimental program was performed at PSI using 
In-house pulsed lasers. 

As a result of the mission euialysis, several attractive missions were 
identified with perhaps the most attractive near term one being the orbital 
repositioning of satellites. In addition, the ed>8orptlon physics scaling ex¬ 
periments served to establish the wavelength scaling of breakdown and absorp¬ 
tion in several gases of interest. For eurgon, the experimental measurements 
were also supported by corresponding theoretical model calculations of the 
plasma formation process. Finally, preliminary single pulse 0.35 vm thruster 
performance data were acquired and compared with numerical code simulations of 
the experiments that model the laser absorption physics, equilibrium gas chem¬ 
istry (including multiple ionization) and gas dynamics expansion. The latter 
data Indicated that laser to blast wave energy conversion efficiencies In ex¬ 
cess of 50% are possible at this wavelength. 

For the present program, work was carried out in three areas; 1) ex¬ 
perimental thruster performance experiments to extend the previous single- 
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pulse 0.35 ym data to multiple pulse operation, 2) modification of the numer¬ 
ical flow code to simulate the fluid dynamics appropriate to repetitively- 
pulsed operation, i.e., a nozzle that has only partly refilled with fresh gas 
prior to each laser pulse, and 3) an analytical assessment of the optical per¬ 
formance of a self-focusing paraboloidal nozzle. 

The thruster performance experiments were carried out using a small 
scale conical nozzle, external focusing optics, emd pulsed CO 2 (X»10.6 ym) and 
pulsed XeF (1*0.35 ym) lasers as the propellant heating sources. All experi¬ 
ments were performed in a vacuum environment. With the XeF laser (Ep*2.7 J, 
Tp*0.6 ys) and argon propellant, a nuucimum specific impulse of 500 seconds was 
achieved with a conversion efficiency of laser energy to propellant )cinetic 
energy upwards of 50%. With the same propellant and the CO 2 TEA laser 
(Ep*>6.5 J, Tp*'3 ys), the maximum specific impulse was 700 seconds with a laser 
to Icinetic energy conversion efficiency ■ 20%. Finally, using hydrogen pro¬ 
pellant and the CO 2 laser as the heating source, the experimental data indi¬ 
cate specific impulse up to 3000 seconds with corresponding energy conversion 
efficiencies upwards of ■■ 70%. The experimental results are presented in Sec¬ 
tion 3. 


In support of the thruster performance experiments, the detailed com¬ 
puter model that was developed previously was modified and used to simulate 
the nozzle fluid dynamics in the conical thruster. The model treats non¬ 
steady nozzle flow of a real gas with laser heat addition. Major enphasis in 
the present work was directed toward simulating the fluid dynamics appropriate 
for repetitively-pulsed operation, i.e., the situation in which the nozzle has 
only partly refilled with fresh gas between successive laser pulses and the 
remaining downstream section of the nozzle contains the low pressure tail of 
the previous blast wave. For a given initial propellant mass distribution and 
laser input, the computer program calculates such quantities as the exhaust 
gas momentum, kinetic energy, and specific impulse. The model calculations 
are described in Section 2. 

Finally, an assessment has been made of optical design criteria for a 
self-focusing paraboloidal nozzle. The primary conclusions reached are: 

1} the mission requirements for off-axis performance allow the surface to be 
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of only very modest optical quality, 2) for off-axis angles I 10°, only very 
modest magnifications (S 25) can be anticipated, 3) optimum optical and fluid 
dynamic performance appear to be mutually exclusive, consequently a compromise 
is required, and 4) the Impact of off-axis focusing and brea)cdovm on thruster 
performance should be evaluated. The details of the optical assessment tas)t 
are presented in Section 4. 
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2. MODELING OF A PULSED LASER-HEATED THRUSTER - SINGLE PULSE 
AND MULTIPLE PULSE OPERATION 

2.1 Introduction 

The main purpose of the modeling effort was to calculate the flow in a 
nozzle which was only partly filled with fresh gas. This simulates the case in 
which the time between laser pulses is short enough so that the nozzle does not 
refill completely after the flow from the stagnation chamber restarts. The re¬ 
gion of the nozzle near the throat contains the fresh gas which has entered 
through the throat. However, the rest of the nozzle contains heated gas left 
over from the previous pulse, which is at much lower density. 

An estimate of the time necessary to fill the nozzle can be obtained 
from the ratio of the nozzle length to the limiting velocity of the incoming 
gas. For 300 K 2 urgon the limiting velocity is 5.59 x 10^ cm/s, so for a 10 cm 
nozzle length it takes 180 us to fill the nozzle. The time for the laser¬ 
generated blast wave to traverse 10 cm might be another 10 us. So if the inter¬ 
pulse time is less than 200 ps (more than 5000 pulses per second) the tall of 
the previous pulse will still be in a 10 cm nozzle when the next pulse occurs. 
Such conditions were explored in the experiments (see Section 3). 

The purpose of having a pulse repetition rate higher than the rate 
which permits complete refilling is to increase the average specific impulse. 

It is clear that the fluid mechanics will be different if part of the nozzle is 
filled with the lower density gas representative of the tail of the previous 
pulse, instead of the fresh gas which has come in through the throat. 

An early study of the fluid mechanics of pulsed laser propulsion by 
Simons and Pirrl (Ref. 2.1) dealt with the multi-pulse problem. However, the 
valuable analytical modelling in that paper was restricted to a perfect gas, and 
utilized blast wave theory, in which the shock wave is always hypersonic 
relative to the gas in the nozzle. In the work described in Ref. 2.2 these two 
restrictions were removed. A computer model of the non-steady nozzle flow of a 
real gas with laser heat addition was constructed, and used to perform single 
pulse calculations. 
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During th« course of the present work, this computer progrem, described 
in Ref. 2.2, was modified to simulate the flow in the partially-filled noczle 
and some runs were made to compare the filled and partially filled cases. This 
part of the work is described in the next section. 

The program was also used to aid the experimental diagnostics by pro¬ 
viding shock wave position vs. time (trajectories) as well as shock peak pres¬ 
sure vs. shock position, for some early experiments in quiescent air. The cor¬ 
relations resulting from these calculations are presented in the second section 
below. 


2.2 Second Pulse Simulation 

Some additional output of the LSDKS computer program was developed in 
order to aid the interpretation of the results of the calculations. To motivate 
the additional output, the quasi-one-dimensional conservation equations which 
are solved are restated here: 


Mass: 


|PA + . 0 

3t 9x 


Momentum: 


dpuA . 9 . 2«. dA 

+ SZ tA(p+pu )] - P ^ 




dx 


Energy: 


1^ [pA(e+uV2)] + 1^ [puA(h+uV2)] 


The equations are written in conservation form. The quantities in the time de¬ 
rivative terms are usually called the conserved quantities, while those in the 
space derivative terms are the fluxes of mass, momentum and total enthalpy. 


The computer program was modified to print and plot the x-dlstributions 
of pA, puA, pA(e-t-u^/2}, and the x-distributions of the time integrals of puA, 
pAu^, and puA(h-tu^/2). The latter three flux integrals might be called the 
mass, momentum and energy fluences. 


The average exhaust velocity for a pulsed propulsion device can be de¬ 
fined in tenns of a mass-averaged exhaust velocity and the mass fluence as 


u 

avg 


2 

/ pAu dt// pAudt 
o o 


( 2 . 1 ) 








where the integrals are evaluated at the exhaust plane. The denominator is the ex¬ 
pelled mass, the numerator the exhaust momentum. 

To simulate the presence of the tall of the previous pulse in the 
nozzle, a reduced density region is required. n>e computer program initially 
fills the nozzle with steady supersonic flow from given stagnation conditions, 
assuming a perfect gas. A simple way to introduce a reduced density region into 
this flow is to merely divide the density by some given factor in the region 
downstream of the desired nozzle station. However, it is not desirable to in¬ 
troduce a pressure discontinuity, since the calculation scheme will treat that 
as an interface to be resolved by another shock wave, which is not the intent. 

To avoid a pressure discontinuity wherever the density is reduced, the internal 
energy can be increased by the same factor (equivalent for a perfect gas to in¬ 
creasing the temperature). This will keep the variation of pressure continu¬ 
ous. 


This idea was Implemented in the computer program by introducing as new 
input quantities a station at which the density and internal energy change, and 
the factor by which the former decreases and the latter increases. 

A number of runs were made in real argon, using the thermodynamic table 
for triply-ionized argon described in Ref. 2.2 which is valid for densities 
above 10”^ g/cm^ and internal energies up to 2.5 x 10^2 erg/g. Below 10"^ g/cm^ 
the program uses the perfect gas law for argon. 

The nozzle geometry was the conical shape, with a 10*^ half angle and a 
throat diameter of 3.1016 cm. 

"nie location of the density and internal energy changes was taken at 2 cm 
from the nozzle throat. For a stagnation temperature of 300 K, which yields a 
limiting velocity of 5.59 x 10^ cm/s, 2 cm is traversed in 36 ps. This condition 
was one chosen for the second pulse experiments, and was used in the calculations to 
simulate those experiments. 

The stagnation conditions used were 10 atm and 300 K. The runs were made 
in the blast wave mode of the program. A given amount of laser energy is put into a 
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small region of the nozzle near the throat to determine the Initial thermodynamic 
conditions there, and then the calculation proceeds to resolve the resulting discon¬ 
tinuity, much like a shock tube flow. For these runs, 1 J was put Into the region 
extending from the throat to 0.1 cm do%mstream. This Is a volume of 1.125 x io~^ 
cm^. The density there was taken to be uniform at the value In the perfect gas noz¬ 
zle flow at 0.1 cm, which Is 3.66 x 10'^ g/ca?. The resulting Internal energy for 
1 J Is 2.429 X 10^^ erg/g, leading to a temperature of 73,000 K and a pressure of 
2.0 X 10® dynes/cm^ (1950 atm). 

The density reduction factor was taken as unity (no reduction), 100 and 
1000 In three runs, to compare the effects of having the tall of the previous pulse 
In the nozzle. In order to see which density reduction factor Is most realistic, 
the density at 2 cm from the throat was plotted vs. time for each of these three 
runs, as shown In Fig. 2.1. The density at that location when the nozzle Is filled 
with gas from the stagnation chamber Is 8.48 x 10“® g/cm®, shown by the line at the 
top. The open squares are for the run with no density reduction, the open circles 
for a factor of 100 reduction and the filled circles a factor of 1000 reduction. 
Figure A shows that there Is not much difference In the density history for the re¬ 
ductions of 100 and 1000. They ]3oth reach 8.5 x 10"^ g/cm® at aJbout 21 ps, and ex¬ 
trapolation Indicates they would reach 8.5 x 10"® g/cm® at about 40 ps. For the no- 
reductlon run, the densities are slightly higher, reaching these values at 24 ps and 
55 ps, respectively. At 36 ps, which It takes to traverse 2 cm at the limiting 
velocity, the density at two centimeters Is nearer to the factor of 1000 reduction 
than the factor of 100, which Indicates that the former Is perhaps a more realistic 
simulation. 

Of course. In the actual case, the properties of the gas In the tall have 
a distribution which Is different from that obtained by merely changing the nozzle 
flow state by some factor. The present calculations should l>e looked upon as a 
qualitative Indication of the effect of having the previous pulse tall in the 
nozzle, not a precise numerical simulation. 

In order to evaluate the effect of the density reduction we will compare 
the results of the computer calculations for the no-denslty-reductlon and factor- 
of-1000 reduction cases. 
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Figure 2«2 shows u^^yg x distributions at various times, where x is cen¬ 
timeters downstreun of the throat, u^yg measured in units of 1.8625 x 10^ cn/s, 
and the times of the various distributions are in ws. Figure 2.2a is for no density 
reduction, while Fig. 2.2b is for the factor of 1000 reduction at 2 cm. The reduced 
density case shows considerably higher values. To convert u^yg to specific impulse 
in seconds, one multiplies u^yg by 190. 

Recall from Eg. (2.1) that at any x and t, the value of u^yg is the ratio 
of momentum flux passing that x, integrated up to t, to the s^une integral for the 
mass flux. Thus, u^yg depends on both x and t. Until the shocic arrives at any 
station, the density is fixed and the velocity is approximately the limiting veloc¬ 
ity, which is 5.5875 x 10 ^ cm/s for 300 K. This is, in fact, the value that u^yg 
takes on in front of the shock, as Fig. 2.2a shows. When the shock passes, the gas 
is accelerated, and u^yg increases. 

In the no-density-reduction case, the shock keeps swallowing significant 
amounts of mass as it moves, and this mass must be accelerated. In fact, the mass 
per unit length the shock swallows is nearly constant, since it is the mass flow in 
the nozzle divided by the undisturbed gas velocity, which is nearly the limiting 
speed for most of the nozzle. As a result, the average velocity decreases with x at 
any given time, and also decreases with time at a fixed x, after the shock passes. 
Figure 2.2a shows that at 10 cm, shortly after the shock passes (32.6 (is) the Zgp is 
190 s. At 16 cm and 61.6 |is, the Zgp has dropped to 170 s. 

When there is a density reduction, the shock swallows only one-thousandth 
as much gas per unit length, after the reduction, and so can move faster and accel¬ 
erate the gas to a higher speed. Figure 2.2b shows that at 10 cm Md 13.8 us, the 
shock has passed and the Igp is 722 s. At 17 cm and 21.8 us, the Igp is down to 
646 s. 

The increase in Igp with density reduction is brought on by the consid¬ 
erable reduction in mass fluence. This effect can be seen in the distributions 
of mass fluence presented in Fig. 2.3. The curves for no density reduction are 
shown in Fig. 2.3a, where the top of the ordinate scale is 1.5 x io~^ in units of 
3.14 X 10~3 g. At 10 cm and 32.6 us, 1.82 x 10~4 g passed, while at 16 cm 

and 61.6 us, 1.48 x 10**^ g has passed. The station of maximum mass fluence at a 
given time is some distance behind the shock wave. For a factor of 1000 density 
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Fig. 2.2a Mass averagad velocity distributions at various times. No density 

reduction. Real argon for stagnation conditions of 300 K and 10 atm 
with 1 J deposited between the throat and 0.1 cm. 
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Fig. 2.2b Mass averaged velocity distributions at various times. Density re¬ 
duced by a factor of 1000 downstream of 2 cm. Real ugon for stagna¬ 
tion conditions of 300 K and 10 atm, with 1 J deposited between the 
throat and 0.i cm. 
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Pig. 2.3a Mass fluence distribution at various times. No density reduction. 

Real argon for stagnation conditions of 300 K and 10 atm, with 1 J 
deposited between the throat and 0.1 cm. 
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Moss Fluence (3.1^1 x 10“^ g) 
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Fig. 2.3b Mass fluence distributions at various times. Density reduced by 1000 
at downstream of 2 cm. Real argon for stagnation conditions of 300 K 
and 10 atm, with 1 J deposited between the throat and 0.1 cm. 
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reduction the corresponding curves are shown in Fig. 2.3b where the top of the or*- 
dinate scale is 3 x 10"^ in units of 3.14 x 10~^ g (which is one-fifth the value 
in Fig. 2.3a). Here the shape of the mass fluence curves is much different than 
in Fig. 2.3a reflecting the fact that nearly all the mass started out between x-0 
and x»2. This amount of mass, 8.98 x 10~^ g, has the value 2.86 x 10"^ on the 
ordinate of Fig. 2.3b, and the figure shows that at 23.8 us, when the shock is 
nearly at 20 cm, only 90% of that amount of mass, at most, has passed any 
station. Some of the remaining 10% is still in the initial part of the nozzle, 
while some may represent numerical error, since the time integrations at each 
station were done rather crudely. The amounts of mass which have passed 10 cm 
at 13.8 US, or 17 cm at 21.8 us are lost at the bottom of the ordinate, but are 
less than 3 x 10~^ g. It is these small values which are the denominator of 
Uavg at the corresponding points in Fig. 2.2b. 

The numerator of Uj^^g, the momentum fluence, is presented in Fig. 2.4. 

For no density reduction. Fig. 2.4a, the top of the ordinate is 0.15 in units of 
5.85 X 10^ dyne-8, while the top of the ordinate in Fig. 2.4b is 0.1 in the same 
units. Again, the values for the no-density-reduction case are higher than 
those for the factor-of-1000 reduction. However, this reduction is not as large 
as for the mass fluence, so the momentum to mass fluence ratio, Uj^yg, is higher 
in the density-reduced case, as we have already seen in Fig. 2.2. 

Figure 2.5 shows that the density-reduced case accelerates the gas to 
higher speeds than the no-reduction case. The gas speed for the no-reduction case 
is presented in Fig. 2.5a where the top of the ordinate scale is 2 in units of 
1.8625 X 10^ cm/s. The rapid oscillations up to 5 cm go with the 61.6 and 62.9 us 
distributions. They are caused by the discontinuity in the gas law at a density 
of 10"^ g/cm^ mentioned above, where real argon is used above 10”^ and perfect ar¬ 
gon below 10"^. Unfortunately, from x»1 to 5 cm the density is very close to 10“^ 
at those times, which produces the large oscillations in the velocity as well as 
all the thermodynamic variables. However, Fig. 2.5a nevertheless shows gas speeds 
after the shock in the range of 0.9 to 1.6. The corresponding curves for the 
density-reduced case are shown in Fig. 2.5b, where the top of the ordinate scale 
is 5 in units of 1.8625 x iqS an/a (which is 2.5 times the value in Fig. 2.5a. 

The gas speeds after the shock are in the range of 3.3 to 4.5, about three times 
as )vigh as for the no-reductlon case* 
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Fig. 2.4b Momentum fluence distributions at various times. Density reduced by 
a factor of 1000 downstream of 2 cm. Real argon for stagnation con¬ 
ditions of 300 K and 10 atm. with 1 J deposited between the throat 
and 0.1 cm. 
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Fig. 2.5a Gas velocity distributions at various times. No density reduction. 

Real argon for stagnation conditions of 300 K and 10 atm, with 1 J 
deposited between the throat auid 0.1 cm. 
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The tick marks on the curves in Fig. 2.5b give the location at which the 
mass, integrated from the throat, equals the mass originally located from x-0 to 
2 cm, namely 8.98 x 10~^ g. This is the gas which represents the new gas in¬ 
jected from the stagnation chamber for the current pulse. The gas from the tick 
marks to the corresponding shock fronts was the reduced-mass gas originally 
downstream of 2 cm, meant to represent the gas remaining from the previous 
pulse. The figure shows that both the residual gas and the newly Injected gas 
have been accelerated to high speeds. 

This simulation shows that if we fill the nozzle only partially with 
fresh gas, 2 md fire the laser while the tail of the gas from the previous pulse 
is still in the nozzle, we can impart a higher momentum per unit mass to the gas, 
and thus a higher specific impulse. 

As an ed.d to the experimental diagnostics, plots were also made of the 
pressure vs. time at two stations in the nozzle where pressure transducers were 
installed, namely 4.4 cm and 9.4 cm. These plots are shown in Fig. 2.6, where 
Fig. 2.6a is for the no-reduction case, and Fig. 2.6b for the factor of 1000 
density reduction. In the former case, we see that the pressure jumps to its 
highest value when the shock 2 urrives, and then decays, for both stations. For 
the reduction case. Fig. 2.6b, the pressure behavior is somewhat different. The 
pressure at 4.4 cm jumps when the shock arrives, but after a small wiggle, con¬ 
tinues to rise for some time before decaying. The pressure at 9.4 cm also jumps 
upon shock arrival, then decays, but reaches a minimum and rises again. (The 
dashed part of the curve is an interpolation necessitated by an error in specify¬ 
ing the location of the pressure station during part of the calculation, so the 
depth of the reduction is not known. However, the right end is defined by two 
correctly calculated points, so there is no doubt the curve rises.) 

The behavior of the pressure shown in Fig. 2.6b has some Implications for 
the interpretation of the pressure transducer data obtained in the double pulse 
experiments (see, for example. Fig. 3.3). In those experiments, the pressure 
peak that is observed at each transducer following the second laser pulse is 
taken to indicate the arrival of the slug of propellant gas heated lay that pulse. 
The propellant exhaust velocity is then Inferred from the transit time of this 
"wave" between the two pressure stations. It appears from Fig. 2.6b that this 
interpretation may have to be examined more closely. 
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Pig. 2.6b Pressure history at pressure transducer locations. Density reduced 
by a factor of 1000 downstream of 2 cm. Real argon for stagnation 
conditions of 300 K and 10 atm, with 1 J deposited between the throat 
and 0.1 cm. 
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2.3 





Calculations for Blast Waves in Quiescent Air 


As an aid to the interpretation of some experiments in quiescent air, 
three calculations were performed. In all cases, the initial state of the gas 
in the nozzle was 1 atm, 300 K, with zero velocity. Energy was deposited into 
the zone near the conical nozzle throat. The gas there was taken to be at the 
same density as the nozzle gas, 1.1719 x lo~^ g/cm^, but at 4000 K. The size of 
the energy deposition zone was varied to allow different amounts of energy to be 
deposited. In one run, the region from the throat to 3 cm was taken as hot, 
corresponding to a deposition of 4 J. In a second run, the region to 1.8 cm was 
hot, corresponding to a deposition of 1 J. In both cases, perfect air, Y*1.4, 
was used. 


The information desired from these runs was the trajectory of the shock 
wave, and its peak pressure. Correlation formulas for these quantities were ob¬ 
tained in a manner similar to that described in Section 3.14 of Ref. 2.2, using 

the appropriate similarity theory of blast waves into quiescent air. The length 
2 1/3 

scale is (Z/Pq&q), where E is the energy deposited, and and ao are the den¬ 
sity and sound speed in the quiescent air. (The latter two for the present case 
are 1.1719 x lO"^ g/cm^ and 3.479 x lo^ can/s.) The time scale is this length 
divided by a^. Therefore, the similarity variables are 

5 - x/(E/p a T - t/{E/p (2.2) 

o o o o 


for length and time. 

The shock trajectory in these variables was calculated for each run. 

The location of the shock at each time was defined by the maximum pressure point, 
and the corresponding value of C was denoted as Cp. A plot of Cp vs. t for the 

4 J and 1 J cases are shown as the circles and triangles in Fig. 2.7. It can be 

seen that the scaling represented by Eq. (2.2) reduces the two trajectories to 
the same curve, in spite of the factor of 4 difference in energy. 

In the experiments to be described later, probes were at x~4 and 9 cm, 

which translate to a pair of values of C for each of the two energies. These 
values are denoted by horizontal ticks on Fig. 2.7. Only the values of C between 
these Buurks are interesting in using the data from the two probes. 
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To correlate these points« a straight line in the log-log plot of 
Fig. 2.7 was fitted to the points by least squares/ excluding the three lowest 
circle points and the one lowest triangle. This portion of the plot is in the 
region where the calculation is just starting, at early times, and shows some 
deviation from a straight line. The equation of the line is 


Cp - 2.76 tO*55926 


(2.3) 


When we revert to physical variables, 
E - 7.16 X 10"^^ X 


using Eg. (2.2), 

6.807^^3.807 

P 


and solve for E, we find 
(2.4) 


in cgs units. This relation expresses the energy in terms of the time of shock 
arrival at the station Xp, and can be used to infer the energy from the measure¬ 
ment of pressure peak arrival time. 


A similar relation laetween the peak pressure wd the location of this 
peak can be obtained. For the pressure, blast wave theory suggests that the 
similarity variable is 


p p o o o 

where Po is the pressure in the quiescent gas, which is 1 atm in the present 
case. A plot of lip vs. C for the 4 and 1 J runs is presented as the circles and 
triangles in Fig. 2.8. Again the two runs fall on a single curve except at low 
values of Cp, which correspond to early times. The vertical ticks again are the 
locations of the two probes for the two runs. 

A straigi.c line was fitted to the circles and triangles, deleting the 
three circles and two triangles at the upper end. The fit is 


TT ^ f p —1.918 

n - 6.6103 C 
P P 


vdilch translates to 


E 


2.55 X 10^ 


•( 2 . 6 ) 


(2.7) 
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in cgs units. This relation can be used to infer the energy from a measurenent 
of peak shock pressure. 

A third run into quiescent air was made as a check case« to compare with 
a run made by Brode in the early days of the use of computers for fluid mechanics 
calculations, Ref. 2.3. He calculated the blast wave resulting from the explosion 
of a sphere of air (Y>1.4) at ambient density but a pressure of 121 atm, into an 
ambient gas at 1 atm. He used Lagrangian coordinates and a calculation algorithm 
introducing artificial viscosity, in contrast to the Eulerian coordinates and Leuc 
calculation algorithm used here. 

This case has an energy deposition which corresponds to approximately 9 J 
in the 10° nozzle used for the present calculations. The results of this calcula- 
tlon, in similarity variables, are shown as the squares on Figs. 2.7 and 2.8. On 
the former, the squares deviate only slightly from the curve defined by the cir¬ 
cles and triangles. In Fig. 2.8, the deviation is larger. However, it must be 
remembered that the initial pressure in the 4 and 1 J runs was only 13.3 atm, com¬ 
pared to 121 in the Brode run. It also appears from Fig. 2.8 that the squares are 
merging with the circles and triangles as the pressure drops. 

The results of the present calculation for the Brode case were compared 
with a small-scale plot of peak pressure vs. distance in Ref. 2.3. The agreement 
seemed to be satisfactory, lending additional confidence to the results of the 
present computer program. 
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3. EXPERIMENTS 


This section describes the experimental measurements that were made to 
assess the performance of a pulsed laser-heated thruster. For the most part, 
the experiments that were carried out were very similar to those performed 
previously^*^ - both in measurement technique emd method of zmalysis. The ex¬ 
periments were performed using a small scale conical rocket nozzle and either 
argon or hydrogen propellant. The propellant gases were heated by pulses of a 
TEA CO 2 laser (X > 10.6 pm, Ep ^ 10 J) and/or an e-beam pumped XeF laser (X " 
0.35 pm, Ep ■> 4 J). The laser radiation was focussed into the nozzle using 
external focussing mirrors and concentrated to sufficient intensity to achieve 
propellant breakdown just downstream of the throat. The resulting dynamics of 
the laser-generated propellant shock waves were monitored using high speed 
piezoelectric pressure transducers mounted in the nozzle wall. Using the nu¬ 
merical code simulations described in the previous section, the pressure 
transducer data were analyzed to infer estimates of rocket performance. The 
performance parameters evaluated include conversion efficiency of laser energy 
to blast wave energy, laser energy to exhaust kinetic energy, and specific im¬ 
pulse, i.e., mass averaged exhaust velocity. All experiments were performed 
with the nozzle in a vacuum background (P < 10*^ atm). 

The following subsections describe the experimental apparatus and 
technique, the rocket performance data and its analysis, and, finally, the 
conclusions that can be drawn from the experimental results. 

3.1 Experimental Apparatus 

A schematic diagram of the overall experimental apparatus is shown in 
Fig. 3.1. In the following subsections, each major component of the apparatus 
will be described. 

3.1.1 Lasers 

The Lumonics K-103 pulsed CO 2 TEA laser was installed in a carefully 
designed and electrically shielded metal box in order to minimize EMI noise 
pickup by the electronic diagnostics. Attention was paid .< shielding against 
rf radiation, avoiding ground loop problems, and filtering the power line 
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Inputs* The CO 2 laser was placed as close to the laser propulsion test chamber 
as space allowed, thus minimtzln? the spreading of. the beam as it travelled 
along the optical train and obviating the need for any beam reducing optics* 

The details of the pulsed XeF laser (0*35 um) have been described pre¬ 
viously*^*^ 


3*1*2 Test Chamber and Optical Train 

The vacuum test chamber used was the seune as that described previous¬ 
ly*^*^ As shown in Fig. 3*1, a KCl window (2" diameter and 1/8" thickness) 
was placed on the plexi-glass flange to allow transmission of the CO 2 laser 
pulse (10*6 um) into the vacuum test chamber. The CO 2 laser beam was then 
turned 90° by a two-by-two inch gold surface mirror and directed to an 8" 
aluminum coated concave surface mirror (f * 160 cm). The two-by-two inch 
folding mirror was placed In such a way that It fit Inside the hole of the 
square annular XeF laser beam* 

Ihe XeF laser was directed, as shown in Fig. 3*1, to the 8" mirror de¬ 
scribed in the above. The two laser beams were then directed onto the 3" 
aluminum coated off-axis parabola* This off-axis parabola mirror focused the 
two laser beams approximately 21 can from the mirror* The focal location of 
the OO 2 laser pulse was carefully identified by observing breakdown sparks in 
air* For the XeF laser a He-Ne laser beam precisely adjiisted to match the 
profile and propagation axis of the XeF beam, was used to determine the focal 
point by observing the He-Ne image* It was found that the focal plane of the 
CO 2 laser beam was approximately 2*5 cm beyond that of the XeF laser* Ibis 
was the result of the finite wavefront curvature of the CO 2 beam* Neverthe¬ 
less, when the apex of the rocket cone was placed 0*6 cm beyond the focal 
point of the XeF laser, the CO 2 laser was powerful enough to induce breakdown* 
As a result, this configuration was used throughout the experiments* * 

To assess how much optical energy was entering in the rocket, it was 
necessary to evaluate the reflectance and transmittance of the various mirrors 
and windows in the optical train* For the optical train associated with the 
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XeF laser beam, the following values were used: quartz window - 0.92 trans¬ 
mittance and 0.08 reflectance: aluminum coated surface mirrors - 0.88 reflec¬ 
tance. The resulting transmission was 0.655. For the CO 2 laser optical 
train, the following values were used: gold mirror, R « 0.97; KCl window, 

T * 0.92; aluminum coated surface mirrors, R > 0.95. Resulting transmission 
was 0.78. 


3.1.3 Rocket Nozzle Design and Propellant Feed System 

A schematic diagram of the rocket exhaust cone and the nozzle Is shown 
In Fig. 3.2. The rocket was made of aluminum 2 md the Inside surfaces were 
polished using conventional polishing and buffing techniques. The throat di¬ 
ameter was 1 mm. Piezoelectric pressure transducers mounted flush with the 
Inside wall surface were at 4.3 cm and at 9.3 cm from the apex of the cone, 
respectively. 

For most of the nuis, argon gas was used as the propellant and gas 
flow was controlled by a fast acting electrically activated solenoid valve. 

The solenoid valve maintained gas flow for approximately 25 ms after it was 
activated, thus minimizing the pressure rise In the time vacuum chamber. A 
piezo pressure transducer was also mounted on the rocket plenum to monitor the 
pressure history as well as the establishment of steady state gas flow through 
the rocket nozzle. This piezo pressure transducer was calibrated against a 
conventional gas regulator pressure gauge. The calibration was accomplished 
by Installing the same solenoid valve and plenum assembly (without nozzle) 
just after the regulator, and then pulsing the gas Into the plenum. In this 
way, pressure drop between the regulator and plenum was minimal. For the ac¬ 
tual -ocket experiments, in which 5 feet of 6 ran i.d. tubing was required be¬ 
tween the regulator and solenoid valve, the maximum plenum pressure was found 
to be approximately 90% of the regulator delivery pressure. 

3.2 Experimental Diagnostics 

3.2.1 Piezo Pressure Transducer 

Two Klstler pressure transducers were used to monitor the laser- 
induced shock wave eurrlval time to calculate shock wave parameters and the 
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post shock pressure history, nie piezo pressure transducers were mounted In 
Teflon® plugs to minimize mechanical vibrations, and they were mounted flush 
with the Inside surface of the conical rocket wall. The two piezo pressure 
transducers vere Installed at 4.3 cm and at 9.3 cm from the apex of the cone. 
The response time was nominally 1 Ps and the responslvlty was 25 mV/psl. 

3.2.2 I.aser Energy Measurement 

The XeF laser energy output was monitored by a Lumonlcs 500 calori¬ 
meter measuring the reflected beam from the quartz entrance window (total re¬ 
flectivity of 8%) as shown In Fig. 3.1. 

Energy measurements were carried out always simultaneously with the 
laser-induced shock transit time measurements, since the energy ouput of the 
XeF laser varied from run to run. This allowed the evaluation of blast wave 
energy or kinetic energy conversion efficiency of the XeF laser energy. 

CO 2 TEA laser (10.6 ym) energy output was monitored using the same 
calorimeter. In a separate measurement from the transit time measurement, 
since the CO 2 laser energy output was very consistent from run to run. 

The CO 2 laser beam was directed to the calorimeter through an 8.75% 
transmittance beam attenuator. The observed energy output was found to be 
consistently 11.9 J. 

3.3 Experimental Technique and Procedure 

The test chamber was evacuated to approximately 2 x 10 '^ torr. For 
most of the experiments, argon gas was used as the propellant and was pulsed 
by a fast solenoid valve to establish steady state flow conditions In the 
rocket nozzle. 17 ms was required to reach the steady state condition over 
the delivery pressure range of from 50 pslg to 150 pslg. The CO 2 TEA laser 
was then pulsed to create a breakdown In the argon gas and to expand the gas 
outside of the rocket cone. The breakdown near the apex was of sufficient 
pressiire to temporarily stop the argon gas flow through the nozzle. After a 
predetermined time delay, the XeF laser was fired Into the rocket. The pre- 













determined time between laser pulses was chosen to be short enough so that 
propellant gas refilled only a small portion of the nozzle. As the XeF laser 
was fired, a simultaneous measurement of energy was carried out as described 
in the previous section. Pressure wave transit times, and pressure history in 
the plenxim and in the rocket cone were recorded by photographing oscilloscope 
traces of the transducer outputs. 

In the double pulse experiments, the purpose of the first laser pulse 
is simply to generate a nozzle flow field for the second pulse that simulates 
what would be encountered in ein ^ thruster. As described above, the break¬ 
down due to the first pulse creates an overpressure downstre£un of the throat 
which temporarily stops the flow from the stagnation chamber. Then, as the 
blast wave expands down the nozzle and the pressure decays sufficiently, the 
flow restarts and the nozzle refills. The extent to which the nozzle refills 
depends on the delay time to the svtbseguent laser pulse. Thus, at the time 
the second (or laser pulse enters the rocket, the region of the nozzle 

near the throat contains fresh, cold gas which has entered through the throat. 
However, the rest of the nozzle contains heated, low density gas left over 
from the previous pulse. It is this flow field situation which is modeled in 
Sec. 2. In the subsequent data analysis, experimental values of specific im¬ 
pulse and laser to propellant kinetic energy conversion efficiencies are de¬ 
rived from second pulse pressure transducer data. 

Other experiments were carried out to determine the conversion effici¬ 
ency of laser energy to blast wave energy. For these, only one laser was 
fired, either the XeF or the CO 2 laser, and the same measurement procedures 
were performed. 

The experimental test matrix is given in Table 3.1. The time between 
laser pulses was varied from 30 to 100 ps. For the shortest time interval, 
the shock due to the deposition of CO 2 laser energy was still in the nozzle. 
The longest time delay was limited by signal to noise and the lower efficiency 
of converting deposited energy to directed kinetic energy. The time Interval 
between laser pulses was controlled a California Avionics digital time de¬ 
lay generator. 
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TABLE 3.1 


Test Matrix of Laser Propulsion Experiments 



Footnote: Numbers in the test matrix indicate mass of propellant (in grains) 

heated by the XeF laser pulse. The propellant used for the above ex¬ 
periment was argon. 
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4 Experimental Results and Data Reduction 
3.4.1 Data Reduction and Error Analysis 

Representative oscilloscope records for the measurements of the XeF 
laser energy, pressure in the plenum, and shocks in the nozzle are shown in 
Figs. 3.3a, b, and c. 

Peak-voltage shown in Fig. 3.3a is proportional to the XeF energy out¬ 
put reflected from the quartz window. The sensitivity of Uimonics SOD calo¬ 
rimeter is 0.39 volt/J and 20% of the XeF laser beam fell outside of the calo¬ 
rimeter, thus the XeF energy output is given by 

E(J) - 12.5 X X 1.25 , (3.4.1) 

0 • 39 

where V is the observed peak-voltage shown in Fig. 3.3. 

Pressure history in the plenum shown in Fig. 3.3b indicates that after 
the solenoid valve opened, the pressure reaches steady state; a large spike 
which appeared at 17 ms shows the time that the CO 2 and the XeF lasers were 
fired. The steady state pressure in the plenum determines the argon mass flow 
rate through the 1 mm dieuneter throat, which is given by 


« * • 

p U A 


(3.4.2) 


- P(atm)(0.2387) 

where p* is the propellant density in the throat, u* is the sonic velocity of 
argon, and A* is the nozzle throat eurea. Mass introduced into the nozzle be¬ 
tween the laser pulses, dm is given by 


Am > m * At , (3.4.3) 

where At is the time between the two laser pulses. Equation (3.4.3) 'is an 
upper limit to the expelled mass since it assumes that the propellant flow 
restarts to the steady state value Immediately after the OO 2 laser was fired. 
Calculations carried out using the flow simulation code described in Sec. 2.2 
suggest that, in fact, for the conditions of the experiment, flow from the 
plenum may not restart for approximately 3 microseconds. Thus, for laser in¬ 
terpulse times ■ 30 microseconds, the error inherent in using Eq. (3.4.3) is 
less than 10%. 
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Fig. 3.3 Represer'-ative oscilloscope 
records of the XeF laser energy 
measurement (-a), the pressure 
history in the plenum (-b), and 
the laser-induced shock wave meas¬ 
urement (-c). (a) 50 ms/div, 50 

mV/div: (b) 10 ms/div 100 mV/div; 

discontinuity on the plateau re¬ 
gion indicates the time when the 
C02 and the XeF lasers were fired: 
(c) 10 ps/div, upper 500 mV, lower 
100 mV/div, respectively; f indi¬ 
cates the CO 2 laser firing; * and 
indicate the arrival of the CO 2 
laser-induced blast wave at the 
first and second pressure trans¬ 
ducer station, respectively; '(j' 
indicates the XeF laser firing; 

A and 5 indicate the arrival of 
the XeF laser-induced shock wave 
at the first and second pressure 
transducer station, respectively. 


'aooinv j 

" • -A 

I'T 

» . 

•• : 


• 

rzn 


•• • 


TTT7 

* 








. 






t.; .. 

— 

^.. 

• ' > * 

1 11 1'' 

itV 

fjlim 

lW4 

r 

r 



1 

Wfd 

Kav 

t ^ 



n 





i: • ■' 

•1 . 


1* 








r.'- 



y 





1 

^ j 


• • r * 

* ' 



R 





/ ! 

* 

m 

V 

□ 





c 






















































































In Fig. 3.3c, the following events were identified as shown with signs: 
(1) the time the CX )2 laser was fired (t): (2) the arrival of the CO 2 laser- 
induced blast wave at the first station (*) and at the second station (•); 

(3) the time the XeF laser was fired (30 ys after (1))('^); (4) the arrival of 
the XeF laser-induced shock wave at the first (A) emd second stations (&)• 

Ihe CO 2 (or XeF) laser energy deposited in the steady state flow of 
argon is given by 


_ 0.4142 * P(atm) * {X(cin)>^*^® .. . .. 

E(J) * -r-rr- (3.4.4) 

{t(us)}^*^® 

where P(atm) is the pressure in the plenum, x(cm) is the distance between the 
breakdown point and the piezo pressure transducer stations, and t(ys) is the 
transit time from the end of the breakdown to the surrival of the blast wave at 
each pressure transducer station. The derivation of formula 3.4.4 is in Ref. 
3.1. Also, for experiments carried out in stationary air, the laser energy de¬ 
posited in the rocket cone is given in terms of pressure by 

E(erg} » 2.41 x lO"^ x pi.5144 first piezo pressure (3.4.5) 

transducer station 


or 

E(erg) “ 2.83 x 10“^ x pi.5144 the second piezo pressure (3.4.6) 

transducer station 

where P is the observed pressure jump from the pre-shock condition (in dyne 
cm*'2). The derivation of these equations is given in Sec. 2.2. 

The laser-induced propellant exhaust velocity, averaged over the two 
piezo pressure transducer stations, w^s obtained by measuring a transit time 
between the two pressure spikes as shown with A and & in Fig. 3.3c. This anal¬ 
ysis was only applied to the second laser heating pulse. 

Similarly the exhaust kinetic energy of the propellant was taken to be 
K.E. - ^ Am • v^ (3.4.7) 
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where v is the average propellant exhaust velocity (as described above) and 
Am is the mass expelled between the first and second laser pulses. Specific 
impulse is taken to be V/g where g is the gravitational constant. 

Since the events of interest in the rocket cone take place on a time 
scale comparable to the response time of the piezo pressure transducers, and 
over distances comparable to the transducer dieuneters and uncertainties in the 
breakdown location, it is important to estimate the maximum errors incurred in 
the experiments. A sample error analysis of the data in Fig. 3.3c is given be¬ 
low. 

Using Eq. (3.4.4) the maximum error in determining the CO 2 laser- 
induced blast wave energy is given by 


0.4142 * {■ 


4.38 Px' 
. 3.38 


3.38 Px 


=* 0.4142 X {2.12 + 3.17 + 0.74} (3.4.8) 

where P = 10.1 atm, t = 16 ms, and x = 9.3 cm and Ax, At, and Ap are taken to 
be 0.3 cm, 1.0 ms, and 0.5 atm. Then, the absolute value of AE is given by 
2.5 J, which corresponds to a i 40% error. As seen in Eq. (3.4.8), the ac¬ 
curacy of measuring the laser energy deposited in the propellant was limited 
primarily by the response time of the piezo pressure transducer and the uncer¬ 
tainty of the precise breakdown location in the rocket cone. Similarly the 
maximum error in the determination of the propellant kinetic energy resulting 
from the second laser pulse is given by 
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where x » 5.08 cm, P * 10.1 atm, t « 10.5 ps, and ■ 30 ps and 

> 1 PS, Ax > 0.3 cm, AP > 0.5 atm, and At > 1.5 ps. Then, the absolute value 
of AE is given by 0.43 J, which corresponds to - 50%. As seen in Eg. (3.4.9), 
the biggest errors arise from the finite dimension of the piezo pressure trans¬ 
ducer head and uncertainties in the transit time measurements. As is illus¬ 
trated above, the systematic errors associated with the present measurements 
are significant, i.e., - 40% for the blast wave energy, t 50 % for the exhaust 
kinetic energy, and f 20% for the specific impulse (average propellant exhaust 
velocity). To reduce these errors, improved diagnostic techniques are needed 
that have improved spatial and temporal resolution. 

3.4.2 Data for Energy Conversion of Laser Energy to 
Blast Wave Energy 

Typical oscilloscope records of the blast wave measurement in station¬ 
ary 1 atm air using the CO 2 and the XeF lasers are shown in Figs. 3.4a and b, 
respectively. For these experiments, the conversion of laser energy to blast 
wave energy was obtained by measuring the elapsed time from breakdown to the 
arrival of the blast wave at each station and measuring the pressure jump at 
each station. Equations (3.4.5), and (3.4.6) were used for the derivation of 
the blast wave energy. The experimental results are shown in Table 3.2. 

Figures 3.5a and b show typical oscilloscope traces used to determine 
the conversion efficiency of XeF laser energy into blast wave energy in an ar¬ 
gon flow. In these experiments, only the transit time from breakdown to the 
arrival of the blast wave at the second station was employed for the calcula¬ 
tion of the conversion efficiency. The conversion efficiencies achieved using 
the CO 2 and the XeF lasers are plotted as a function of mass flow rate in 
Figs. 3.6 and 3.7. It is seen that for the range of mass flow rates investi¬ 
gated the coupling of the CO 2 laser (10.6 wm) beam with argon gas is.relative¬ 
ly independent of the flow rate while the coupling of the XeF laser (0.35 um) 
beam appears to increase with flow rate. Nevertheless, for mass flow rates " 
1.0 g/s, both lasers appear to be converted to propellant blast wave energy 
with an efficiency ■ 50 ± 20%. 
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TABLE 3.2 


Comparison of Two Different Types of Blast Wave Energy Conversion 
Measurements Using the CO 2 e-nd the XeF liasers (in Stationary Air) 


X = 4.3 cm 
X * 9.3 cm 


Observed Energy (J) and Conversion 
From Transit Time 

4.41 J 48% 

3.42 J 37% 


Efficiency (%) of the CO? Laser 
From Pressure Jump 
6.49 J 70% 

5.19 J 56% 


X ^ 4.3 cm 
X = 9.3 cm 


Observed Energy (J) and Conversion Efficiency (%) of the XeF Laser 
From Transit Time From Pressure Jump 

0.677 J 3i% 0.974 J 55% 

0.832 J 47% 1.27 J 72% 
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3.4.3 Data for Conversion Efficiency of Laser Energy to Gas 
Kinetic Energy Using Two Laser Pulses - Determination 
of Specific Impulse 

Typical oscilloscope records of the measurements used to determine 
propellant kinetic energy were shown in Figs. 3.3a, b and c. The propellant 
kinetic energy appearing in the exhaust is calculated using Eq. 3.4.7 with the 
propellant exhaust velocity determined from transit time measurements of the 
propellant pressure wave induced by the second laser pulse (XeF laser). The 
reduced data are plotted in Figs. 3.8-3.12 as the conversion efficiency of 
laser energy to propellant kinetic energy vs. propellant mass. The propellant 
mass plotted is the mass expelled between the first and second laser pulses, 
i.e., m A t, where At is the laser interpulse time. For these experiments, 
the XeF laser energy entering the rocket was typically 2.6 J. 

In Fig. 3.13, we plot the average conversion efficiency achieved vs. 
argon mass flow rate. For each mass flow rate, the conversion efficiency 
plotted represents an average for the various laser interpulse times studied. 
As can be seen, for mass flow rates i 1.6 g/s, the conversion efficiency is 
relatively independent of mass flow rate and appears to take on an average 
value of 25 to 30%. 

From Figs. 3.8 through 3.12 it can be seen that specific impulse of up 
to 500 s was achieved in argon propellants heated by the pulsed XeF (0.35 urn) 
laser. Furthermore, this specific impulse was obtained with conversion 
efficiencies of laser energy to propellant kinetic energy of up to 30 to 50%. 
These results compare favorably with similar results that were presented in 
Ref. 3.1 for argon propellants processed by two successive 10.6 um laser 
pulses. In the latter studies, maximum specific impulses of 500 s in argon 
and 1000 s in hydrogen were achieved with conversion efficiencies of laser 
energy to exhaust kinetic energy of up to 40 to 50%. 

Since in both the previous studies^*^ and the present work, the maxi¬ 
mum specific impulse achieved in argon was 500 seconds, the question arose as 
to whether this represents a fundamental physical limit in argon. However, 
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ment, which was evaluated using Eq. (3.4.9) 
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because the earlier experiments were performed with a nozzle geometry signifi¬ 
cantly different from the present one, i.e., a rapidly expanding paraboloid 
vs. the present 10* cone, it was difficult to draw any strong conclusions. 
Thus, it was decided to perform two pulse experiments in the conical nozzle in 
which both the first and second pulses were CO 2 laser pulses. In principle 
this would allow direct comparison of 0.35 urn and 10.6 urn performance under 
the same nozzle conditions, as well as comparisons of the performance of the 
parabolic and conical nozzles under the same conditions. 

In Fig. 3.14 we plot the specific impulse data obtained for two 10.6 
pm laser pulses heating argon propellant in the conical nozzle. As was de¬ 
scribed above, the specific impulse is determined from transit time measure¬ 
ments of the nozzle pressure waves induced by the second laser pulse. The 
data of Fig. 3.14 are plotted as the specific impulse observed vs. Am, the 
propellant mass heated by the second CO 2 laser pulse. As before, we have 
taken Am to be equal to ^t, where At is the time delay used between laser 
pulses. The diagonal lines plotted with the data represent loci for different 
values of exhaust kinetic energy, i.e., « 1/2 AmV^ « 1/2 Am(glgp)2. As can 

be seen, the maximum specific impulse achieved in these experiments was ap¬ 
proximately 700 s. For a propellant mass of ■ 5 x 10~^ grams, this value of 

specific impulse indicates a maximum exhaust kinetic energy of " 1.2 J. With 

the total CX )2 laser energy entering the rocket being approximately 6.5 J, this 
corresponds to a kinetic energy conversion efficiency of 18.5%. It is inter¬ 
esting to note that a specific Impulse of 700 s in argon is 40% higher than 
any of the previous results. Comparing conical nozzle results only, this im¬ 
proved performance is probably due to the higher laser energy available at 
10.6 lim vs. 0.35 pm, i.e., 6.5 J vs. 2.7 J. Presumably similar results could 

be obtained with an XeF laser if the laser energy were increased. 

Based upon the most simple arguments, the specific Impulse is expected 
to scale as ~ *^Eiaser* that this is not quite the case when we compare 

the specific Impulse obtained using the CO 2 laser and that achieved with the 
XeF laser. From the ratio of the available laser energies, ~ 6.5/2.7, the 
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’resent Conical Nozzle Data 
d* = 0.1 cm 

• m = 2.^1 g/s 
41 m = 1.63 g/s 
g| ffi = 1.37 g/s 

A m = 0.99 g/s 

T m = 0.70 g/s 


Parabola Data (Ref 3.1) 
d* = 0.2 cm 

O rn = 1.1 g/s 
^ ni - 4.8 g/s 


.V,—-5-5 5 j— 

10'5 10 ^ 

PROPELLANT MASS (dm) - grams 

Fig. 3.14 Specific impulse data obtained for argon propellant heated by two CO 2 
TEA laser pulses. Specific impulse is determined from transit time 
measurements of pressure wave induced by 2nd laser pulse. Dashed 
lines plotted with data indicate loci of constant propellant kinetic 
energy (h Am(g Isp)^)- ^ is the mass expelled through throat 

between first and second laser pulses. 











above argument would lead us to predict a specific impulse Improvement of 55%. 
The observed improvement was only 40%. It is clear that in future studies it 
would be worthwhile to investigate how the achievable specific Impulse varies 
with input laser energy. 

It is also interesting to note the higher specific Impulse achieved 
with the present conical nozzle compared to the previous parabolic nozzle. We 
believe the improved performance results from conditions in the conical nozzle 
being more conducive for efficient laser energy deposition. Apparently, at 
the lower mass flow rates, the very rapid initial divergence of the pareibollc 
nozzle leads to downstream gas densities that are too low to support strong 
LSD wave cibsorption. 

Finally, measurements similar to those presented for argon in Fig. 

3.14 were also obtained for hydrogen propellants. The results are given in 
Fig. 3.15. Note the significant increase in specific Impulse that occurs when 
hydrogen is used, l.e., up to 3000 seconds in hydrogen vs. a maximum of 700 
seconds in argon. Qualitatively, this result is not surprising in light of 
the much lower molecular weight of hydrogen. In fact, for specific impulse 
scaling as the inverse square root of molecular weight (assuming constant 
stagnation temperature), an increase of /40/2 • 4.5 is expected. The observed 
improvement is a factor of ~ 4.3. 

Comparing the results in Figs. 3.14 and 3.15 we also note that the 
conversion efficiency of laser energy to propellant kinetic energy is consid¬ 
erably higher in hydrogen than in argon. In fact, at maximum specific im¬ 
pulse, the conversion efficiency achieved in hydrogen appears to be approxi¬ 
mately four times greater than in argon, i.e., ~ 74% vs. ~ 19% in argon. 
Because the data in Fig. 3.6 Indicate CO 2 laser energy deposition efficiencies 
in argon of 40 to 50%, it would appear that the low overall conversion to ex¬ 
haust kinetic energy in this propellant is limited primarily by loss mech¬ 
anisms other than imperfect absorption. Energy loss channels that are likely 
to be important are plasma radiation losses and losses to residual enthalpy of 
the exhaust gas. 
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Present Results^ Conical Nozz 
A m = 0.26 g/s 

• rli = 0 . 5 ^ g/s 


□ Previous Results (Re 
Parabolic Nozzle 
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PROPELLANT MASS (Am) - grams 

Fig. 3.15 Specific inpulse data obtained for hydrogen propellant heated by two CO 2 
TEA laser pulses. Specific impulse is determined from transit time meas 
sureinents of pressure wave induced by second laser pulse. Dashed lines 
are loci of constant propellant kinetic energy (*j Am(g Isp)^)- Am is 
the mass expelled through throat between first and second la^ser pUl^es. 














As was pointed out with the argon results, the present results for hy¬ 
drogen using the conical thruster exhibit a significant improvenent in perfor¬ 
mance over the previous results obtained with the parabolic nozzle,In the 
latter experiments, the best performance obtained with hydrogen propellant was 
a specific impulse of 1000 seconds at a conversion efficiency of laser to ex¬ 
haust kinetic energy of " 40%. We believe the better performance achieved in 
the present investigations, i.e., 3000 seconds with conversion efficiencies 
to ■> 70%, is primarily the result of more efficient laser energy absorption in 
the conical thtirster. In the previous parabolic thruster, the initial nozzle 
divergence was so large that attempts to increase specific impulse hy reducing 
the amount of propellant mass heated led to nozzle gas densities that were too 
low to support efficient laser absorption. 


3.5 Summary of Results and Conclusions 

The results of the present thruster performance experiments can be 
summarized as follows: 

(1) Using a conical nozzle of 10* half angle and external focusing 
optics the following rocket perform 2 mce parameters have been achieved: 

• Laser to blast wave energy conversion efficiencies in eurgon propellant 
of ~ 40 to 60%. Such energy deposition efficiencies were achieved 
using either a CO 2 TEA laser (1 “ 10.6 ym, Tp ■ 3 us, Ep “ 6.5 J) or a 
xenon fluoride exclmer laser (X » 0.35 um, Tp ■ 0.6 us, Ep ■ 2.7 J). 

• Specific impulse in argon of 500 s using the XeF laser (Ep » 2.7 J) 
and 700 s using the CO 2 laser (Ep " 6.5 J). These specific Impulses 
were achieved with corresponding conversion efflciences of laser ener¬ 
gy to exhaust kinetic energy of up to 40% ± 10% and 19% ±5%, respec¬ 
tively. 

• A maximum specific Impulse in hydrogen propellant (using the CO 2 laser 
as the heating source) of ~ 3000 seconds. This specific Impulse was 
achieved with a conversion efficiency of laser to exhaust kinetic 
energy of up to ~ 70%. 

For all the above experiments, specific impulse was determined from trMslt 
time measurements of the pressure wave induced by the second laser pulse in a 
two pulse heating sequence. Typical interpulse times were 20 to 100 us. 


60 








(2) Performance obtained with the conical nozzle and OO 2 TEA laser 
pulses represents a significant improvement over similar experiments performed 
previously^* ^ using a parabolic nozzle. Previously the maximum s^^cific im¬ 
pulse that could be obtained was 500 s in argon and 1000 s in hydrogen. This 
is to be compared to the present results of 700 s and 3000 s, respectively. 

The improvement is believed to result from better laser energy absorption that 
occurs in the conical nozzle when operating at low mass flow rates. 

Based upon the experimental results described above, the following 
conclusions can be made: 

(1) There appears to be no fundamental reason why pulsed laser-heated 
thrusters c<m not be operated efficiently at vlsible/UV laser wavelengths. 

The present results indicate conversion efficiencies of laser energy to pro¬ 
pellant energy at X « 0.35 imi that are comparable to those obtained for X > 
10.6 pm. The only significant question that remains with respect to wave¬ 
length scaling is whether, if a self-focusing nozzle is to be used, sufficient 
focusing can be achieved at the short wavelengths to achieve rapid laser- 
induced breakdown. Optical design considerations for a self-focusing nozzle 
are discussed in Sec. 4. 

(2) It appears that the values of specific impulse that %»ere obtained 
previously using a parabolic nozzle,^*^ i.e., 500 seconds in argon and 1000 
seconds in hydrogen, were limited by poor laser energy deposition at low pro¬ 
pellant masses. Apparently, the very rapid Initial divergence of that nozzle 
created downstream gas densities that were too low to efficiently absorb the 
incident laser radiation. In contrast, the present conical nozzle, with its 
10* half angle divergence, eliminates this problem by maintaining higher gas 
density downstream of the throat. The most draunatic improvement is seen in 
the case of hydrogen for which the achievable specific inqpulse Increased from 
1000 seconds to some 3000 seconds. To achieve similarly high specific lBg>ulse 
in the parabolic nozzle will probably require operating at higher mass flow 
rates and with significantly larger laser pulse energies. Thus, for full 
scale operation, parabolic nozzles %rlll probably not present a limitation to 
rocket performance. 














(3) Finally, It is worth noting the similarity between the rocket per¬ 
formance parameters we have obtained for a pulsed laser-heated thruster and 
those typical of an arc-heated rocket engine. According to Refs. 3.2 and 3.3, 
arc-heated rocket engines typically operate with a specific impulse in the 
range of 400 to 2000 seconds and a maximum propellant temperature of ~ 5500*C. 
The typical working fluid is hydrogen. The internal efficiency of such en¬ 
gines, i.e., fraction of electrical input energy converted to kinetic energy 
of ejected mass, is typically 15 to 20 percent, (but can probably be increased 
to some 50 to 60 percent if the propellant is used as a regenerative coolant). 
The similar performance observed for a pulsed laser-heated thruster is not 
surprising in light of the similarity of the propellant heating process, i.e., 
gas breakdown and plasma formation followed by a thermodynamic nozzle ex¬ 
pansion. 


3.6 Recommendations for Future Studies 

Perhaps the most significant result of the present studies is the ex¬ 
perimental evidence for achieving a specific impulse of some 2500 to 3000 sec¬ 
onds with hydrogen propellant. Unfortunately, because this result depends on 
a theoretical interpretation of the laser-induced pressure waves measured in 
the nozzle, its value is in question. For verification, additional experi¬ 
ments are required in which a more direct, reliable measurement of specific 
impulse is carried out. Ideally, a diagnostic technique is needed that can 
provide a measurement of exhaust gas velocity. 

One approach which appears feasible (particularly for exhaust gas vel¬ 
ocities in excess of 10^ cm/s, i.e., Xgp > 1000 seconds) is to measure the 
Doppler shift of spectral line emissions emanating from the exhaust. In hy¬ 
drogen, a good spectral line to monitor might be the well-known red line, Ha, 
of the atomic hydrogen Balmer spectrum. For particle velocities of 3 x 10^ 
cm/s (Igp “ 3000 seconds), the Doppler shift would be ^p/c or 1 part in 10^. 
High resolution spectrometers are readily available which can resolve 1 part 
in 105. 


More reliable measurements of rocket performance would also be made 
possible by performing longer duration laser tests. Sufficient test time 
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could enable the use of standard rocket diagnostic instrumentation such as 
load cells, thrust stands, strain gauges, ballistic pendula, etc. Of course, 
to obtain specific impulse from measurements of thrust and total Impulse re¬ 
quires that the propellant mass flow rate and total expelled mass be known ac¬ 
curately. 


Finally, independent of performing any new experiments, the existing 
experimental data should be analyzed in greater detail through use of the com¬ 
puter simulation code described in Section 2. Using the multiple pulse code 
with time dependent laser absorption and real gas thermodynamics, calculations 
should be carried out for the repetitively-pulsed performance of the conical 
nozzle operating with argon and hydrogen propellants heated by laser pulse en¬ 
ergies and with laser interpulse times that simulate the experiments. It will 
also be extremely valuable to incorporate into the computer code a subroutine 
that calculates plasma radiation losses. 
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4. OPTICAL DESIGN CONSIDERATIONS FOR A SELF-FOCUSING NOZZLE 


We consider in this section the optical performance of a paraboloidal 
self-focusing nozzle of utility for laser propulsion. The primary conclusions 
are: 

* The requirement for off-axis performance allows the surface 
to be of only very modest optical quality. 

• For off-axis angles • lO^’, only modest magnifications (£ 25} 
can be emtlcipated. 

* Optimum optical and fluid dyn 2 unic performance appear to be 
mutually exclusive, consequently a compromise is required. 

• The Impact of off-axis focusing and breakdown on thruster 
performance and surface erosion should be evaluated. 

Details of the analyses are given in the following paragraphs. 

4.1 Aberration Characterization 

Two general classes of aberrations can be distinguished; determinis¬ 
tic and random. The former are those defects in the image caused by the de¬ 
parture of the actual wavefront from a spherical wave for a mathematically 
perfect reflective or refractive surface. The latter are caused by departures 
of the actual surface from the desired mathematical surface which are random 
because they are associated with manufacturing tolerances. For emalysis pur¬ 
poses. it is convenient to consider each type separately because of the dif¬ 
ference in their basic character. 

4.1.1 Random Aberrations 

Consider the simple optical system shown in Fig. 4.1 which consists of 
a lens of aperture D and focal length Z. Ihe image of a point object at in¬ 
finity is given by 


Kx) 


(Xz)"^ / dvw(v)U(v)exp[- X 



(4.1) 


where X is wavelength. U(y) is the complex wave in the apert\ire plane and W(]r) 
is a simple aperture function (-1 inside and 0 outside the aperture). We now 
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assume that U(v) is characterised by the product o£ a deterministic portion 
(constant for a perfect aperture) and a random phase distortion (- exp i^(v))* 
The average image plane intensity can be shown to be equal to 

x*A 

<I(x)> - P(-)Ij,(x) + —^ [i-pO] / dATjj(A)p^(A)e ~ ~ ( 4 . 2 ) 

(^s) 

where <..•> represents a statistical average, Iq is the deterministic image, 
and E is the energy collected by the aperture. Tq is the optical transfer 
function of the deterministic image (normalized Fourier Transform of 1 q) and 
Pq is the correlation function of the random ccxnponent of the field. 


p(A) - 

<exp i(^(v) - ♦(v~A)]> 



p(A)-p(“) 


Pn(A) - 

[ i-p(«) ^ 

(4.3) 


1 1^ 1 1^ 


P(») - 

|<exp i4>(v)>| 5 |*(1) , 



where $ (1) is the first order characteristic function of the random vaurlable 
♦(v). 


In many cases of interest, the correlation function is much narrower 
than the deterministic transfer function leading to the approximate result. 

)oc 

<I(x)> - {l*(1)l^ P(x) +7 [l-|*(1)|^]s (—)} , (4.4) 

(Xz)^ ~ * * 

when P(x) is the normalized point spread function of the deterministic image 
and S is the spectrum of the random field (Fourier Transform of Pij). 

From this result, we note that energy is conserved (l.e., inttegral of 
<I(x)> over X is equal to E). Furthermore, the energy in each term is given 

by 


(Ed/E) - [♦(Dl^ 
(Er/E) - 1-|*(1)|2 
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(4.5) 









where Eq and are the energies of the deterministic (first term of Eg. 4.4) 
and randm (second term of Eg. 4.4} components of the field, respectively. 

The on-axis intensity is given by: 


<I(o)> -{14(1)1^ + [i-|*(1)|^](A/d)^} , (4.6) 



These general characteristics of the image are shown in Fig. 4.2. 

These results indicate that, in principle, a deterministic image is always 
present (of size Xz/D if unaberrated). However, its energy is attenuated by 
the presence of the random component of the field. The random component of 
the field results in an average intensity which appears as though it were 
generated by a lens of aperture A, the correlation scale of the random field. 

The characteristic function depends on the type of statistics involved 
and the rms phase distortion (o^). Values for several common cases are shown 
in Table 4.1. As can be seen, one wave (2ir) of random aberration substantial¬ 
ly eliminates the deterministic component in all cases. For normal statis¬ 
tics, a guarter and (1/20) wave of aberration results in a 92% and 9% reduc¬ 
tion in energy, respectively. Conseguently, high guality optical systems are 
usually specified at (X/20) or better. 

4.1.2 Deterministic Aberrations 

Except for diffraction, the on-axis performance of a mathematical 
paraboloidal mirror is perfect based on third order aberration theory. For 
off-axis cases (i.e., direction of object not parallel to axis of rotation) 
Coma and Astigmatism occur (Spherical Aberration is always zero). The angular 
diameters of the blur circles associated with these effects cure given by: 


Coma: 



(4.7) 
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TABLE 4.1 First-order characteristic functions for 
various distributions. 


Distribution 

Characteristic function 



Value 

2v 

at a 

* 4ir 

normal 

2 

exp(-o.) 

10-'" 

10-®^ 

uniform* 

2,,, 1/2 , 
sine [(3) a ] 

9 

8 X 10'^ 

10-^ 


4, .1/2 

-4 

-7 

Simpson 

sine [(3/2) ' 0 ^] 

♦ 

2 X 10 

2 X 10 


2 1 

-2 

-3 

exponential 

(1-a^)-1 

V 

2 X 10 

6 X 10 

Laplace 

(1-0^/2)-2 

9 

2 X 10 ^ 

10-* 

gamma 

-oj 

2 V 

io-'2 

OD 

1 

O 

Rayleigh^ 

^F^[1;1/2;-0^/(4-n)] 

io-'5 

10-"" 


+ [ o^/( 4-iT)exp-'0^/( 4-it) ] 

9 9 




*sincx = (sinx/x)* 

^lF^(a; P; i X) ■ confluent hypergeometric function 


■ I where 

n-0 

(a)ti“(a+l) ••• (o+n-1), (o)q“ 1 . 
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(4.7) 


Astigmatism: B, *■ 

A 

where 6 is the off-axis angle, F is the F-number of the system, z is the dis¬ 
tance from the paraboloid apex to the focal point and L is the distance to the 
aperture stop (see Fig. 4.3). 


For the case of interest here, a physical aperture stop is not lilcely 
to be used. Consequently L is defined by the physical extent of the mirror or 
the beam size. In this case, the blur circle di^uneter due to Astigmatism is 
given by 


2F 2 '•' * 

16F 


(4.8) 


The expressions given above are reliable for modest apertures or fields, how¬ 
ever, they break down at larger values. A very conservative estimate of the 
resultant blur due to both effects is given by the sum. We will use a root 
sum square which is perhaps more realistic. 


Inspection of the above forms lead to the obvious conclusion that for 
off-axis performance, a larger F-number is desirable. Paraboloids of various 
F-number and fixed diameter are shown in Fig. 4.4. As the F/Ho. increases 
with a fixed entrance diameter, the paraboloid becomes shallower emd the focal 
point moves farther away from the apex. Above some critical value, the 
resulting paraboloid will no longer function adequately as a rocket nozzle. 

The gas flow will separate from the walls and/or the gas expansion rate will 
not be appropriate. Furthermore, for values in excess of F/0.25, the focal 
point will lie outside the entrance aperture. Unless some type of mechanical 
extension is used, this will likely yield an inefficient motor because the 
thrust will not be vectored within a small enough angle. Consequently, nozzle 
design should be based on a trade-off between the conflicting requirements for 
optical vs. fluid dynamic performance. 


The paral?olic rocket used in laboratory experiments (Fig. 4.5) cor¬ 
responds to approximately F/0.05. Ray traces for this thruster at 0 > 0° and 
15° are given in Figs. 4.6 and 4.7, respectively. As can Ise seen, the beam 
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4.6 Ray trace for F/0.05 paraboloid with 0 























does not focus to any great degree at 6 * 15°> Ray traces for an F/0.1 systeo 
are shown in Fig. 4.8. 


An approximate upper limit to off-axis performance can be obtained by 
the ad-hoc criterion that the aberrated bean (rhs of Eg. 4.7 and 4.8 converted 
to spatial size) must be smaller than the physical dimension of the parabola, 
i .e., 


16F 




(4.9) 


The limiting angle defined by this expression is given in Fig. 4.9 as a func¬ 
tion of F-number. The magnification associated with this spot size is also 
given in Fig. 4.9. 


The results of Fig. 4.9 indicate that for low F-number systems (< 0.1), 
the laser beam will not focus for off-axis angles greater than a few degrees. 
Optical performance at larger angles can only be determined by a detailed ray 
trace which considers multiple reflections or via experiment. Restricting the 
system to small angles has significant implications in terms of overall mission 
configurations. 


4.2 Discussion 

The analysis given in the prior section indicates that because of 
fluid-dynamic considerations, a low F/No paraboloid will be required. This 
will result in relatively large off-axis deterministic aberrations which will 
produce a large focal spot size at any angle above a very modest value. This, 
in turn, reduces the magnification dramatically. Consequently, any detailed 
system study should include the requirement that the energy density of the 
input laser beam be within roughly a factor of (10-25) of that required for 
propellant breakdown. 


When coupled with a maximum off-axis angle of 10*^ or less, these two 
requirements inould seem to be rather restrictive in terms of mission configura¬ 
tions. Obviously, these requirements can be relaxed implementing a larger 
F-number system. As mentioned previously, this has an effect on perforaance as 
a t)iruster. 
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Fig. 4.9 Limiting off axis emgle and corresponding magnification for various 
F/# paraboloids. 
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Magnification a 0 . 








To establish a complete map of allowed operational pareuneters will re¬ 
quire a determination of the limiting F-number Imposed by the fluid dynamics. 

While the leurge spot size is^lled by off-axis performance limits the 
achievable magnification, there is one advantage in terms of optical quality. 
Referring to Fig. 4.3, we see that the random image size is given approximately 
by (Xz/A) where A is the correlation scale of the surface. Provided this 
dimension is smaller than the spot size associated with the deterministic off- 
axis performance, there is no reason to require a high optical surface quality 
corresponding to a substantial fraction of energy in the diffraction image 
(X/10 or better). Data from one vendor of paraboloidal mirrors suggests that 
this scale is of order 0.1 mm which is consistent with other data relative to 
surface (non-optical) finishes. For a one meter, F/0.05 thruster, this value 
yields a random image size of less than 1 mm at a wavelength of 0.35 y which is 
much smaller than the 20 cm nozzle size at nominal focus (for a one meter aper¬ 
ture). Consequently, several waves of random aberration could he tolerated 
without serious impact on performance. To determine the actual amount allowed 
would require an alternate analysis which takes into account surface slope 
(i.e., o^/A), and shadowing effects which can become important for low inci¬ 
dence angles on a surface with steep slopes. The analysis of Section 4.1 does 
not adequately account for these effects, i.e., it is appropriate for a moder¬ 
ately good surface. This alternate analysis is available in several standeurd 
texts but it is probably not necessary to pursue this issue in detail. 

One effect which has not been evaluated is the impact of off-axis 
focusing and breakdown on thruster performance. As the beam moves off-axis, 
the focal point also moves off-axis. This may be important. However, as in¬ 
dicated in the prior results, the beam cpiickly becomes large and so it may not 
be appropriate to conceptualize the situation in terms of breakdown occurring 
at a small point which is located at a substantial distance from the eucls. 

A related issue is surface erosion which may occur as the focus and 
breakdown region moves close to the walls. Again, however, it is probably in¬ 
appropriate to consider the problem in this fashion. Rather a more conserva¬ 
tive approach would be to determine the effect on the surface when the break- 
do*m occurs in a large volume which has a dimension of order of the nozzle in 
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the vicinity of nominal focus. It seems reasonable that this assumption viould 
lead to a reasonable estimate of surface erosion effects. 

In summary, it appears that the allowed off-axis angle is significant¬ 
ly limited the low F-numbers which probably are required for suitable fluid 
dynamic performance. This has the further effect of yielding large spot sizes 
which limit the achievable magnification. This does have benefit in that only 
very modest optical surface quality is required. To establish a complete map 
of allowed operational conditions will require the evaluation of the maximum 
allowed F-number which is controlled by the fluid mechanics. Off-axis focus¬ 
ing could be a problem, and a conservative approach to evaluating surface 
erosion would be to assume that breakdown occurs in a volume with dimension of 
order of the nozzle size in the vicinity of nominal focus. It should be noted 
that these conclusions are based on a rather simple analysis. Depending on the 
input beam profile, adequate energy densities which will lead to gas breakdown 
may exist in a small region even for large angles and nominal spot sizes. This 
effect could be significant and hence, a detailed ray trace which includes beam 
profiles should be Implemented to provide a complete assessment of these 
issues. 




